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Abstract 

We propose a matrix evolution equation in (x, A;)-space for flavour singlet, unintegrated 
quark and gluon densities, which generalizes DGLAP and BFKL equations in the relevant 
limits. The matrix evolution kernel is constructed so as to satisfy renormalization group 
constraints in both the ordered and antiordered regions of exchanged momenta k, and 
incorporates the known NLO anomalous dimensions in the MS scheme as well as the NLx 
BFKL kernel. We provide a hard Pomeron exponent and effective eigenvalue functions 
that include the nj-dependence, and give also the matrix of resummed DGLAP splitting 
functions. The results connect smoothly with those of the single-channel approach. The 
novel Pqa splitting functions show resummation effects delayed down to x = 10"'^, while 
both Pga entries show a shallow dip around x = 10~^, similarly to the gg single-channel re- 
sults. We remark that the matrix formulation poses further constraints on the consistency 
of a BFKL framework with the MS scheme, which are satisfied at NLO, but marginally 
violated by small n^/A'^^-suppressed terms at NNLO. 



-"^On sabbatical leave of absence from Dipartimento di Fisica, Universita di Firenze and INFN, Sezione di 
Firenze. 



1 Introduction 



Small-x QCD evolution, historically based on DGLAP [1] and BFKL [2] dynamics, has been 
widely investigated in the past years [3-9], leading to a better understanding of the two ap- 
proaches just mentioned and to robust resummed predictions for the gluon density and splitting 
function [10-17]. There is now a remarkable consensus [18] among various resummation ap- 
proaches on the resulting gluon evolution kernel, and a satisfactory comparison of some of them 
to experimental data [15, 19]. 

The basic idea underlying the progress of the resummation approaches just mentioned, 
lies in the observation [10] that the BFKL kernel embodies an infinite number of subleading 
contributions which are collinear singular. These terms are parametrically large and need to be 
taken into account in order to achieve consistency with the rcnormalization group (RG). The 
techniques for incorporating such terms differ in detail according to the various authors, but 
lead eventually to similar results. 

However, all approaches developed so far limit themselves to a consistent resummation 
scheme only for the evolution of the gluon density. The quark-sea contribution is obtained on 
the basis of fc-factorization of the qq dipolc [5] in the DIS factorization scheme and/or by adding 
the qq contribution to the next-to-leading-log x (NLx) BFKL kernel. The main purpose of the 
present paper is to devise a resummed small-x evolution scheme in a coupled matrix form, so 
as to treat gluons and quarks on the same footing and in a collinear factorization scheme which 
is as close as possible to a predetermined one, e.g. the MS scheme. 

Our matrix approach, in the collinear limit, has the advantage that it complies automatically 
with the matrix factorization of the integrated partonic densities in the singlet evolution, and 
thus is able to incorporate the known low-order anomalous dimensions for any value of a; = 
N — 1, the moment index. On the other hand, in the high-energy limit, the (gauge- invariant) 
unintcgrated partonic densities arc well defined by fc-factorization around different values of uj. 
To be precise, the gluon unintcgrated density is defined around uj = and the quark around 
a; = — 1. Therefore, in the leading high energy region — that is around a; = — we are able to 
include the known LLx+NLx BFKL kernel in the gluon channel only, thus leaving the quark 
entries somewhat unconstrained from the fe- factorization standpoint. 

Note however that assuming a BFKL framework in matrix form is a demanding requirement, 
because fc-factorization implies resummation formulae for the anomalous dimension matrix up 
to NLx level. Therefore, incorporating both exact low-order anomalous dimensions (say, in the 
MS scheme) and NLx expressions (from the exact BFKL kernel) imposes on our matrix kernel 
some nontrivial consistency relations expressing the requirement that collinear and high-energy 
schemes do not conflict with each other. They are discussed in detail in the following, and 
we flnd that they are satisfied in the MS scheme up to NLO level, while they are marginally 
violated at NNLO, by small terms in the gq entry of relative order Uf/N^. For this reason we 
restrict ourselves, in this paper, to the NLO-NLx level. This means that, starting at NNLO, 
our splitting functions will be in some matrix scheme which is not the MS scheme. Higher order 
effects on the scheme change have been studied in [14] and have been found to be comparable 
to rcnormalization scale uncertainties. 

Besides the NLO-NLa; information mentioned above, we impose the general requirement of 
consistency with the renormalization group in both ordered and antiordered configurations of 
exchanged partonic momenta fc's. This is best expressed in the 7 <-> l+cu — j symmetry of the 
kernel, where 7 is conjugated to logfc^. We enforce this symmetry by the so-called consistency 
constraint [20-22] which introduces an a;-dependence in the (leading) kernel, so as to resum 
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those parts [9] of the higher order BFKL kernels which are required by the RG. This procedure 
follows previous papers [10,11] and is used in particular for the gg matrix element of the kernel. 

Despite all such requirements, there is a considerable ambiguity in our approach which 
is tied up to the matrix structure, because the 7 and uj dependences of the various matrix 
elements of the kernel are constrained only to a limited extent by the collinear and high-energy 
limits. We therefore introduce in sec. [2] some further requirements, mostly related to the pole 
structure of the kernel in the 7 and uj variables, by requiring it to have at most simple poles. 
This assumption is quite natural in the case of u; — > because of the BFKL limit, and follows 
by the cu-expansion method [10] in the 7^0 case. The leading twist pole structure of the 
kernel at lowest order in is basically K ~ Fq [(I/7) + 1/(1 + u; — 7)], where Fq denotes the 
LO DGLAP anomalous dimension matrix. The full kernel to second order in is constructed 
in sees. |2] and [3] according to the requirements stated above. It also contains running coupling 
effects at the scales suggested by the NLx BKL kernel and by the RG. 

In the frozen case we calculate in sec. S] the resulting anomalous dimension matrix, 
and its eigenvalues 7-i-(q;s, tu). The leading eigenvalue at high energies, 7+, contains important 
resummation effects in the as/uj variable. We also obtain resummation formulae for Tqg and 
Fpq, the latter being specific to our matrix approach and not directly obtained on the basis 
of the NLx BFKL kernel only. In sec. |5] we present results for the effective eigenvalue (or 
characteristic) functions uj = x±(q^s,7) as inverse functions of the anomalous dimensions, and 
for the hard Pomeron exponent ujs{as). In the case with running coupling described in sec. El we 
provide the resummed DGLAP splitting function matrix in x space, obtained by the numerical 
deconvolution method proposed in [11,23] and generalized to the matrix case. Details of the 
matrix kernel and of the anomalous dimension expressions are left to Appendices A-C. 

2 Basis of matrix formulation 

The general purpose of this paper is to provide integro-differential matrix equations for uninte- 
grated parton distributions, which interpolate between DGLAP evolution equations [1] in the 
hard scale variable log and the high-energy BFKL evolution equation [2] in the rapidity-like 
variable logl/x. Despite the high-energy and collinear factorization constraints, the above 
interpolation is subject to considerable ambiguities, due to the following facts: 

(a) Off-shell, unintegrated densities are defined in a gauge-invariant way by fc-factorization of 
gluon and quark exchanges around different values of the moment index uj = N—1, namely 
a; = for the gluon and u = —1 for the quarkj^ Therefore, in the high-energy region 
— that is around the leading value uj = — only some effective gluon equation (which 
incorporates the high-energy quark contributions) is constrained by the BFKL limit, 
whose kernel has been calculated perturbatively [2,4-7]. This makes the interpolation of 
the kernel to generic uj values more ambiguous for the quark entries. 

(b) The collinear limit constrains the matrix kernel for both quarks and gluons, but in a 
factorization-scheme dependent way, and only in the strongly ordered region of transverse 
momenta ■ ■ ■ ^ ^ k2 ^ ■ ■ ■ and in the anti-ordered one. This limit only restricts 
the singularities of the kernel in the variable 7 (conjugated to logfc^) so as to reproduce 

^This is because of the spin 1/2 of the quark exchange, which leads to an energy dependence of the cross- 
section of type at high energies. For alternative approaches to the definition of unintegrated densities 
see [22,24]. 
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the low order anomalous dimension matrix. In addition, the collinear ^ anti-collinear 
relationship implies the existence ofa7<->-l+u; — 7 symmetry, whose form is however 
quite general, depending again on the factorization-scheme. 

2.1 Basic criteria for the kernel construction 

In order to tame the ambiguities of the off-shell continuation mentioned above, we shall use a 
few basic criteria which — we shall argue — can be consistently imposed and correspond to a 
factorization-scheme choice for both high-energy and collinear limits. 

Let us refer to a matrix kernel Kaf,(Q;s, cu), acting on fc-space, such that the parton Green's 
function is given by 

Qabiuj;k,ko) = [l-K{as,uj)]^j^ {k,ko) , {a,b = q,g). (2.1) 

In the frozen as limit, the kernel matrix elements are diagonalised in 7-space and given by the 
eigenvalue function Kafe(as, cu, 7). Our first basic assumption is that in the collinear limit 7 — >■ 
and uj fixed, the matrix kernel K shows simple poles only, in the form of a 7-expansion 

K = i K W (as, uj) + K(i) (as, u) + j K^^) {a,, + (7') . (2.2) 

Here, the I/7 singularity is natural because of the DGLAP limit, and would be the only term 
present in a pure evolution equation in logfc^. In fact, it implies (sec. 12.21) that the one-loop 
anomalous dimension matrix ro(a;) is given by K^q\uj), the coefficient of the I/7 pole of lowest 
order in as (while the higher order terms r„ involve K^^\ K^'^\ • ■ ■ as well). 

Note however that higher powers of as/7 could have been present alsojf] and do actually 
occur in the normal formulation of the NLx BFKL kernel [6,7]. By eq. (12. 2p we explicitly 
exclude such possibility in our matrix kernel, while the BFKL kernel will be recovered by 
proper algebraic manipulations (sec. [3]). 

Our second assumption is analogous to (12. 2p with uj and 7 interchanged. In the high-energy 
limit of a; ^ with 7 kept fixed we require simple pole singularities in the a;-expansion 

K = - oK(as, 7) + iK(as, 7) + ^ 2K(as, 7) + M , (2.3) 

UJ 

where, in addition, 

oK,, = = oK,g . (2.4) 

The 1/uj singularity is natural because of the BFKL limit and would be the only term present 
in a pure evolution equation in logl/x. It implies that the eigenvalue function of the LLx 
BFKL kernel is given by Xo(7) ~ 0^0(7), the coefficient of the 1/uj singularity of lowest order 
in as (while the NLa; BFKL kernel, discussed in sec. [3l involves iK(7) also). 

Therefore, higher order singularities in a; — which are present in the anomalous dimension 
at higher order — will be obtained (sec. H]) by using the rough anomalous dimension relation 
7 ~ tts/^ in the subleading 7-dependence. The fact that only gK^g and oKgg possess the 1/u 
singularity is related to the fact that in usual factorization schemes [8], only Tgq and Tgg show 
a LLx dependence on the Os/uj variable. 

■^For instance, the rough hhx anomalous dimension relation 7 ~ as/o; could be replaced in the subleading 
(jj-dependence, thus producing higher powers of as/"f. 
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There is a third important assumption, which deals with the relationship between collinear 
and anti-collinear orderings of exchanged transverse momenta. Both orderings are to be in- 
corporated in our off-shell formulation and simple kinematical considerations show that the 
variable conjugated to logfc^ in the reverse ordering is 1 + a; — 7. Therefore, values of 7 and 
1 + — 7 must be related by some symmetry, and we shall assume, out of simplicity, 

Kab{l,Uj) = Kab{l + UJ--f,Uj) . (2.5) 

It is perhaps useful to recall the formal basis for the symmetry ( 12. 5p . Let us write the 
fc-factorization formula for the A + B ^ X differential cross-section in the form [25] 

M; = /2^U) "^'^MSikM'AK,.), (k^lkl), (2.6) 

where we have lumped in the A,B superscripts the dependence on the hard scales Qa,Qb of 
the process. Then the change of energy-scale from kko to, say, can be incorporated by the 
change of kernel 

K^'''\k,k';u;) = (J-^ K['='=«](fc,fc';cj) , (2.7) 

or, at frozen a^, by the cu-dependent shift [7, 9] of the corresponding eigenvalue functions 

X^''\j,u;) = x^''-\j-^,u). (2.8) 

On the other hand, in the one-channel case — namely when one considers only gluon dynamics 
— the A ^ B symmetry of ctab implies the k ^ k' symmetry of the kernel K(A;, k'; u) and the 
7 1 — 7 symmetry of the eigenvalue functions x^'^'^^Kl ^ ^) ■ Therefore, at energy-scale /c^, the 
7 ^ 1 + — 7 symmetry of eq. (12. 5p holds for x^'''^\l:^), whose superscript will be dropped 
from now on. 

In the matrix case, the thorough discussion of sec. 12.31 shows that the collinear <-> anti- 
collinear symmetry of the matrix kernel is expected to have the more general form 

K{l + uj- 7, uj) = 5(cu)K^(7, uj)S-\uj) . (2.9) 

Therefore, eq. (12.51) is obtained by choosing the similarity transformation S so as to have 

SK^S-^ = K , (2.10) 

and represents yet another restriction of our off-shell schemeU 

In the following we shall show in more detail how to construct the matrix kernel so as to 
satisfy the known collinear /high-energy limits with LO-LLx accuracy (sees. [2^ 12.31) and NLO- 
NLx accuracy (sec. E]), within the scheme restrictions provided by assumptions (12.21 12. 3^ 12. 5p . 
We note from the start that eqs. (12. 2p and (12. 3p impose consistency relations on the anomalous 
dimensions, which show up in a novel NLx resummation formula for Tgg, to be discussed in 
detail in sees. 14.21 and 14.31 

^We use the symbol x to denote eigenvalue functions of kernels considered in or related to previous works 
on small- a; resummations in the gluon-channel. We keep the symbol K to denote both kernels and eigenvalue 
functions specifically designed for this matrix formulation. Note also that the x's are perturbative coefficients 
of expansions in ds = asC^/7r, therefore differing in normalization by a factor (2Ca)~'-"^^^ from their K 
counterparts in eq. (|2.1ip . because d = as/27r. 

^The choice of eq. (|2.5p must be supplemented by a corresponding choice of impact factors in order to satisfy 
the symmetry for observable cross-sections. In the realistic NLO-NLa; case the similarity transformation S in 
eq. p.lOp is expected to be an operator in fc-space also. 
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2.2 Form of kernel at LO-LLa? accuracy 

In order to discuss the above features in more detail in the frozen ctg hmit, we introduce the 
triple expansion 



7 



a 



n,m,p=0 



2tt 



(2.11) 



where pKi™'' are 2x2 matrices in the a = q,g indices, and we note that m,p>0, consistently 
with the simple pole assumption of eqs. (12. 2[ 12. 3p . We also use the notation K^'^\as,uj), 
pK(Q!s,7), Ki™'''(a;) and pK„(7) to mean partially resummed coefficients, as already done before. 
In this paper we limit ourselves to two terms (n = 0, 1) in the frozen ag-expansion, which will 
be able to accommodate the LLx and NLx BFKL kernels. However, running coupling effects 
will be introduced by various scale choices for the various terms (sec. 13.31) and this implies in 
general an infinite series when expanding around a fixed scale. 

Let us first show how to construct Kq so as to be consistent with the collinear and high- 
energy limit at LO-LLx accuracy. We denote by F the anomalous dimension matrix, with the 
expansion 



n=0 



n+1 



TJuj) 



(2.12) 



where we recall the small-a; behaviour (Ty = T^nj = nf/2) 



( 



ro(^) 



4T 



2C. 



0(1) 
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^ + (cu) 
+ 0(1) 



(2.13) 



and, in the MS scheme. 



ri(^) 

with the two eigenvalues 




9CfCa - AOTfCp (12C^ - A6CA)Tf 



(2.14) 



7+,o 
7-,o 



2C 



— + (1) , 7+,i 

UJ 

ATfCp 



27> 

9uj 



IOC, 



13 



0(1] 



3C 



+ 0M 



(2.15) 
(2.16) 



A 



We then write the generalised BFKL equation for the unintegrated parton densities J-'i{k] u) 
in the form 

jr = KJ^ + jrsource ^ ^2.17) 

where the source jF^°"'''^<= is local in fc-space. It is then straightforward (cf. sec. 14.11 and app. El) 
to derive, for frozen Og, DGLAP type equations for the integrated densities 



(2.18) 
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of type 



where 



/■ - gliji = E . (2.19) 



To = Ki°^ (2.20a) 
ri = KS°^ + Ki'^KS,°) (2.20b) 
= k(°) + K^C + K^'^Kf) + K?)(0^ + (K«)^C , (2.20c) 

and so on. This identification rehes on the expansion fl2.1ip . which in turn rehes on the assumed 
single 7-pole structure of fl2.2p . see app. \M 

We note that this procedure imphes that the part of the anomalous dimensions, proportional 
to K(^) with 2 > 0, at a given order of perturbation theory can be generated from the lower 
orders. 

Eqs. f l2.20p can be used to constrain recursively the 7^0 singularities of Kq, Ki, ■ ■ ■ for 
a given set of low order anomalous dimensions, for instance in the MS scheme. In particular, 
it determines Kq°^ = Fq, as noticed after eq. (12. 2p . but does not fix Kq^'', which is therefore a 
scheme-changing parameter. At LO-LLx level, we choose the parameterization 

/r,,,o(^)Xc(7) r,,,o(^)xc(7) \ 

Ko(7,^)= , (2.21) 

Vg,,o{^)xn^) ^Xi^(7) + [r,,,o(.;) - ^] xni)J 

where |^ 

X^(7) = 2^(1) - ^(7) - ^(1 + a; - 7) , (2.22) 

reduces to the leading BFKL eigenvalue function in the uj = limit (ipix) is the Digamma 
function), and 

X^Al) = - + T-^ , (2.23) 

7 1 + — 7 

is a simple interpolation of 7 = and 1 + — 7 = poles. Higher twist terms will be needed 
at the NL level in sec. El 

Let us note that in eq. (12.210 we have already incorporated assumptions (12.30 and (12.50 at 
this level. In fact, the form (I2.2ip is consistent with eq. (12.30 because of the BFKL limit of Xo 
in eq. (12.220 and of the absence of 1/u; poles in Kgg and Kqg. Furthermore, the symmetry (12.50 
is present by construction. On the other hand, the simple form of Xc ^q. (12.230 and the fact 
that it appears unchanged in the qq, qg and gq entries and in the part of the gg entry which 
has no l/c(j singularities are all "off-shell" features which we choose for simplicity reasons. 

The fact that Kgg and Kgg are proportional to Xc helps to satisfy the momentum conservation 
sum rule. Indeed, since 

F,,,o(cu = 1) + rg,,o(^ = 1) = Tgg^oiuJ = 1) + Tgg^oiuJ = l) = , (2.24) 

we obtain 

+ = J-, + ^, = X/.t-F,, (2.25) 



^We will use the uj superscript to denote w-shifted [11] eigenvalue functions x'^il) = Xl(7) + X-l(1 + ~ 7)j 
XL being the "left projection" of the eigenvalue function xil) = Xl{"/) + Xl(1 ^ 7) with singularities in the 
half-plane 5R(7) < 1/2 only. 
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where all quantities are evaluated at u; = 1 and Xht is some higher-twist kernel, having sin- 
gularities at 7 = —1, — 2, ■ ■ ■ . It is then easy to show that the sum-rule violation is at most 
(ofg), instead of (a;s)|j This feature will be improved in sec. Oby modifying the parameter- 
ization (12.211) at NLO-NLx level so as to reduce the violation to (a^). 

2.3 General form of the coUinear ^ anti-collinear symmetry 



• • • • • • 



Figure 1: A sequence of splittings going towards smaller xasa—^b— >■€—>■ d, which may be 
collinear, anti-collinear or some mixture. 

The role of the symmetry (12.51) in implementing RG properties deserves a special discussion 
for the matrix kernel. In fact, its typical effect is to produce two pole terms in 7 and 1 + 0; — 7 
(see eq. (12.231) ). which are supposed to describe the correct product of F-matrices for both 
direct and reverse orderings. The situation is pretty clear for coUinearly ordered particles. In 
fact, referring to fig. [1], with x decreasing from right to left we have a sequence of splitting 
functions such as 

...F,,F,,F,,... h <^^<^b<^^ (2.26) 

describing a going to b, b to c and so on, as predicted from K also. On the other hand, in the 
anti-collinear limit the DGLAP splitting functions need to account for the opposite splittings, 
b to a, etc. 

... FerfFb^F,, ... > ^'>> (2.27) 

"'d ^ "'c ^ ^ "'a 



(ndc(rOc5(r^ 



Iba 



Then, one would naively expect that the anti-collinear pole in the kernel be associated with 
F^. This seems to work fine as long as we consider a complete chain of anti-collinear splittings. 
Problems arise however when trying to join collinear and anti-collinear chains. Firstly there is 
an issue of colour factors: the anomalous dimension Tij implicitly includes a factor Ni for the 
number of varieties of parton i that can be produced {Nq — 1 if z is a gluon, 2nfNc for an 
(anti)quark). For each exchanged particle in fig. [T]that factor should be included exactly once. 
In the collinear limit it is included in the branching to the right of a given exchange (e.g. for b 
it is included in F^a), while in the anti-collinear limit, as written in (I2.27p . it is included to the 
left (in The) ■ If we are to consider a single evolution from right to left containing both collinear 
and anti-collinear splittings we should ensure that the Ni factors are consistently included to 
one side, for example in the branching to the right of the exchange. One then needs to correct 
the splitting function (F^)jj for an anti-collinear i —>■ j splitting by a factor Nj/Ni. 



^In the single-channel case (fq — 0) eq. (|2.25p would imply that the sum rule violation is higher-twist only. 
In the matrix case, a higher twist component is expected on top of the perturbative component discussed here. 
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The second issue that arises relates to high-energy factorization. For each exchanged gluon 
we have a factor Xjio. In a sequence of coUinear branchings that factor is associated with the 
splitting function to the right of the gluon exchange (e.g. if 6 is a gluon then it is included in F^a), 
while for anti-collinear branchings it comes from the splitting function to the left (i.e. from V^^. 
This causes problems if we have an anti-collinear splitting to the left of a gluon exchange and 
a coUinear one to the right, since hoth will include a 1/a; factor for the intermediate exchanged 
gluon. However it is necessary for the gluonic part of our Green function to be consistent with 
high-energy factorization, which systematically assigns an exchanged gluon's l/w divergence to 
the larger- a; part of the diagram, i.e. to the right of the exchanged gluon in the coUinear limit. 
Therefore in the case of an anti-collinear branching we should multiply (F^)jj by a factor /i//; 



3 



where fg = 1/uj, so as to ensure that j is never associated with a I/lu factor, while i has it 
when z is a gluon. Note that fq is arbitrary (other than that it should be a non-zero constant 
for cij — > 0) since high-energy factorization is not defined for quarks around uj = — we shall 
discuss its choice below. 

The outcome of this discussion is that colour factor and the high-energy factorization cor- 
rections can be combined by introducing a similarity transformation matrix 

/ 2nfNcfg{uj) \ 

S=\ , (2.28) 

and defining a 'refactored' splitting function matrix F for anti-collinear splittings in an evolution 
that will combine both coUinear and anti-collinear splittings: 



ST'S 



Tc-l 



' '''' Cf fgiuj)'- 91 » 



A matrix kernel will therefore have coUinear and anti-collinear structure of the form 



(2.29) 



K ^ - + , (2.30) 

7 1 + cij — 7 

and will satisfy the coUinear ^ anti-collinear symmetry in the general form (12. 9p . The fact 
that the diagonal entries (in particular the gg element) of F and F are identical is consistent 
with our expectation that the single-channel [n^ = 0) limit should coincide with BFKL, which 
is symmetric in 7 ^ 1 + — 7. The structure of colour factors and Xjuj ensures that chains 
containing coUinear and anti-collinear splittings will have the expected sets of colour factors 
and overall Xjuj factors. 

Note finally that we can further specify S so as to satisfy the symmetry in the form (12. 5p . 
Since 



^gqfl Cp W + 3 ' 

we can simplify eq. (12.301) by exploiting the arbitrariness of fq{uj) and setting it to 



(2.31) 



2T — 
/,(^) = — ^ =^ F = F, (2.32) 



thus providing, at leading level, a fuUy symmetric coUinear structure 

'1 1 

^7 1 + a; — 7 

as assumed in eq. (12.51) . 



K ^ F f- + ) , (2.33) 
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3 The kernel at NLO-NLcc accuracy 



3.1 General structure of Kq 

Let us recall that, while the LO anomalous dimension matrix and the LLx expression of Tgg 
are factorization-scheme independent, the NLO, NLx expressions do depend on the scheme 
(except possibly for the eigenvalue 7+ in the frozen as limit). This opens up the possibihty 
of constructing the kernel so as to reproduce the NLO, NLx anomalous dimensions in a given 
scheme, say MS scheme. However, we have to comply with the restrictions i\2.2\ \2.3\ 12. Sp . in 
particular the requirements of a simple 7-pole structure at fixed u, a simple cu-pole at fixed 
7, and absence of l/a; singularity in Kgg and Kgg. This means, at NLO, that the dg/uj terms 
of Tqq and Tgg cannot be reproduced by an explicit 1/u term in Ki, but should result from 
the 7-dependence of Kgg and Kgg o, where Kq^-*, Kq^"*, ■ ■ ■ are free scheme choice parameters. 
In other words, we have to adjust the sub leading 7-dependence of Kq so as to reproduce the 
known MS anomalous dimensions at NLx level in the form 



(i)n 
Jqg 



for 



UJ 



. 



(3.1) 



The above discussion shows that we have to change the parameterization (12.211) at next-to- 
leading level so as to allow a more general subleading 7-dependence. We choose the following 
one 



Ko(7,u;) 



/rgg,0(^)X^(7) rg,,0(u;)X^f (7) + Agg (uj) X^,{^) \ 

^Tgg,oiuj)x-Al) r,,,o(^)x^^(7) + ^ [Xo (7) - Xc (7)]y 
where Xhtil) a higher- twist kernel possessing the symmetry (12. 5p . e.g. 



(3.2) 



Xhtil) = 3 



+ 



1 + 7 2 + ct; — 7 



xSt(O) = 1 , 



(3.3) 



and Agg is an tu-dependent coefficient which we require to be regular for 3?(a;) > —1 and 
vanishingly small as cu — > 00, e.gjf] 



Agg{uj) = 6gg A{uj) = 6 gg ' 3 



1 



1 + uj 2 + uJ 3 + uj 
'(1) 



A,M 



(3.4) 



The form (13.20 allows one to choose (Kq )gg so as to reproduce the known NLx expres- 
sions (I2.13P and (I2.14p of Tgg and Fg^ in the MS scheme, up to order al/uj. Note that the 
logic is here reversed with respect to the DIS scheme, in which (Ko)gg is directly calculated by 
fe-factorization (see [8]) and NLx resummation formulae for Tgg and Fg^ are derived. By the 
7^0 behaviour 

(3.5) 



7 



Ko - -K^ + K« 



we derive, from eq. (I2.20p for Fi, that = 
In more detail, the NLx coefficient org^.i 

by 



8T//9. 
= lim,„_ 



.0 ^ 



Fgg 1, according to eq. (I2.20bp is given 



Kf) + K^^^K 



(l)l/(0) 



-I qg 



■ I k(°) 4- k(^) k(°)' 



-I qg 



(3.6) 



^This particular choice is motivated by the fact that the z-space function Agg(z) ~ (1 
for z 1 (cf. eq. (|C.6p ). thus not disturbing the large- a; behaviour of the model. 



2)^ rapidly vanishes 
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By expanding the matrix products in terms of their matrix elements and by taking into account 
the conditions (12.41) . only the last term in the r.h.s. of eq. (13. 6p does not vanish, yielding 



or 



.of^o )gg 



[r,,,o(0)c,(0) + 5,gxSt(0)] = 2C, 



3^/ + ^gg 



(3.7) 



where we used the explicit expressions (13. 3p . (13.291) and (l3.3Qbl) . The MS-scheme value or^^,! 
40T/Ca/9 of eq. (12.141) is then recovered provided 



gg 



82/ 
9 



(3.8) 



Note that the corresponding expression for or^g,! involves {si^o'^gg — §f (oKo°'')gg, thus respect- 
ing the colour charge relation T^J^^' = §^r^g^^ > 1) which is apparent in eq. (I2.14p . 

One could repeat the above procedure in other factorisation schemes as well. In the DIS 
scheme, e.g., the value o^qg,i = 52T/C^/9 is recovered by setting 



6. 



DIS 

gg 



UTf 



(3.9) 



3.2 BFKL limit and general structure of Ki 

Having fixed Tgq and Tgg at NLO-NLx level, the remaining constraints (exact NLO anomalous 
dimension matrix and exact BFKL kernel at NLx level) are fixed by a proper choice of Ki. 
Because of the expansion in eq. (12. 3p 



K = - oK(as, 7) + iK(as, 7) + ^ 2K(as, 7) + M 

00 



(3.10) 



where qK has only gg and gq entries, the resolvent 

1 



;i-K) 



-1 



'1 



iK) 



-1 



UJ 



K(l-iK)-^ + 0(cu) 



-1 



(3.11) 



has the gg matrix element proportional to a pure log 1/x evolution form, with the kernel 

[oK(l-iK)-i], 



Tv-BFKL - Tv-BFKL , -2 tv-BFKL 



'99 



(3.12) 



where a= 



asNc 



. We thus arrive at the identification 



(2C^)iro^^^L ^ ^^Ko),, , (2CA)'i^r^^ = (oKi + oKo iKq),, , (3.13) 



which parallels eq. (12.200 for the perturbative expansion, with the difference that it concerns 
the gg entry only, as is appropriate to the fc-factorization of gluon exchange. Eq. (I3.13P is used 
— as in the single-channel case [11] — to derive (oKi)^^ from the known expression of the NLx 
BFKL kernel irf^L and of the Kq kernel Explicitly 



)Ko iKo)gg = 2CaXo(2CaXo + AggXc) + 2CFXc{,TqgfiXc + SqgXht] 



ui=0 



Xo(7) = d^xoil) 



uj=0 



y(l-7) 



(3.14) 
(3.15) 



^Of course, running coupling contributions to Kf^^^ — explicitly considered in sec. 13.31 — are to be sub- 
tracted out. 
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where Agg{uj) = rggfi{u!)—2CA/^^ is the regular part of the gg anomalous dimension. This result 

reduces to the corresponding one-channel [11] nj = subtraction 2C aXoC^C aXo + ^gg^~°^Xc) 
after setting Tgg = Agg = 0. 

We encounter at this point a consistency relation on the factorization scheme for Fi, due to 
the fact that we want to incorporate both Fi and Kf^^^ in a kernel satisfying the simple-pole 
assumptions of sec. 12.11 as better discussed in sec. 14.31 Note in fact that, given Kq, eq (12.201) 
determines K^""* in terms of Fi and eq. f l3.13p determines (oKi)^^ in terms of Kf^^^. Therefore, 
[oK[°'']^^ — the 7-pole part with the l/oo singularity — is determined in two independent ways, 
which should provide the same result, in the form 

o[F,,,i - (K«C),,] = o(7+,i) - (oKo iKo)S^ , (3.16) 

where we have used the fact that the simple-pole part of f^ppKL pj-g^jic^^g ^^j-^g NLx part of 7+,i, 
that is 

o(7+,i) = (2C^)^(i^r^^)^°^ . (3.17) 

We show in sec. ( 14. 3p that the consistency equation (13.16^ is identically satisfied, provided 
oTqg^i ~ [Kq^'*]^^, as assumed in eq. (13. ip . Therefore, both DIS and MS schemes satisfy eq. (13.160 
and can be accommodated by a proper matrix kernel at NL level. 
More precisely, we start from the "Ansatz" 

K.h,.) = r.(.) + (2C,f (i - ^) (° ^,^°^,) . (3.18) 

where 

(i) the function Xi (7) at = is equal to 

~.=0 _ ~ _ oKgg.l _ ^BFKL .0 T a) 

' " {2CaY ^ ^ 

[see eq. ( I3.28P for an explicit expression] and is extrapolated to generic uj values by the uj- 
shift [11] procedure of left and right projections as follows 

Xiil) = [xMl) + [Xi]l(1 + ^ - 7) , (3.20) 
so as to satisfy the symmetry (12.50 : 

(a) in order to minimize momentum sum rule violations, we have added to the high-energy 
pole 1/uj in front of Xi ^ low-energy term —2/(1 + u): their sum vanishes at a; = 1; 
(in) the Fi(a;) matrix is fixed by matching K^^-* + Kq^-^Kq^^ to the known NLO MS splitting 
functions (cf. eq. i^^^)): 

f M = rr'(-) - CC' - (2C.)= (i - (° . (3.21) 

xf*^ being the coefficient of the simple pole at 7 = of Xi (7)- We note that Xc ^l- (I3-I8p is 
again chosen for simplicity reasons according to the symmetry (12.50 . 
The final form of the next-to-leading matrix kernel is then 

K, = (f, - Ki^r) + i2CAr (--T^](' ImJ ■ (3.22) 



u i + uj \o xi-xf^x:. 
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3.3 Running coupling features 

We shall choose the running coupling scales of our kernel as in [11]. We thus associate as(?^) 
[q = k — k') to the LLa; BFKL kernel Xo (t); ^ind Q;s(fc>) (A;> = max{fc, k'}) to all other ones. 
The choice of the intermediate gluon momentum transfer is suggested by the NLx BFKL kernel 
itself, which contains the beta-function dependent term 



"2 (Xo + Xo) 



/ 



6 = 

12 3Ca 



(quoted for a renormalisation scale choice /x^ = k"^), corresponding to the fc-space kernel 



(3.23) 



ir™'^(fc, k') = -b 



(3.24) 



reg 



where the regularization procedure is explained in [11]. Since the term in eq. fl3.24p is accounted 
for by expanding asiq"^) up to NL order, the expression in eq. fl3.23p should be subtracted out 
from the NL kernel considered before. More precisely, the kernel Kf^^^ in eq. fl3.13p is meant 
to have the eigenvalue 



BFKL 
A-l 



Xi - 2^0X0 



+ 2 (^0 + Xo) 



(3.25) 



where the expression in square brackets is the eigenvalue function at energy-scale fc^, obtained 
by the cu-expansion of eq. (12. 8p . and xi is the customary NL eigenvalue [6,7] at energy-scale 
kko, given by the expression 



b 11 
Xiil) = -2[Xo(7) + Xo(7)] - ^Xo(7) - ^ 



TT 



sm 7r7 



cos7r7 
1-27 



3+1 + 



2Tf\ 2 + 37(1-7) 



CI 



'67_7r2 
"'" 36 12 



n=0 



5T 



9C/ 



Xo 7 + 3 + — 

2 4 sm TT'-f 



V^(n+ 1 + 7) - V^(l) V^(n + 2 -7) -^(1) 



(n + 7)^ 



+ 



(n 



if 



It follows that the overall kernel has the structure 



K(fc, k'; u) = a{q^ 



2Ca /O 



,0 



+ a{kl 



+ a^ikl) 



1^(1) k(o) 



ir,^ + {2CAy - 







1+UJ 



;i + 27)(3 -27). 



(3.26) 



(3.27) 



,0 K'^ 



x^K 



Kq, K'^ and K^^^ are the fc-dependent kernels corresponding to the characteristic functions Xo , 
Xc and Xht respectively, and their explicit expressions in (x, /c)-space are provided in app. [C] 
(cf. eqs. f lC.7p . (10.40 and (IC.SP ). The eigenvalue of the gg entry of Ki at = 0, thanks to 



^"its generalisation to include variable renormalisation scale is constructed as follows: single powers of Ofg 

undergo the transformation a{q^) a{x^^q'^) + /?od^(a;^g^) Inx^ and a{ky) a{xf^k'^) + /3oti^(a:^fc?,) Ins^, 

2 ^ 



a{x^^q'^) + /?od^(a;^g^) Inx^ and a(fc?.) 
with /3o — 'iCAb, while quadratic powers of as are modified as d (fc?.) — > a'^{xj^ky). 
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eq. (13.251) and to the subtraction procedure in eq. fl3.14p . is provided by the expression 

Xi 



Xiil) + \ [Xo(7) + Xo(7)] + ^^o(^)— 2^ 



Agg{id)l2CA Cf/Ca 



7(1-7) 7(1-7) 



qgfi 



(0)/2C, 



+ 



7(1-7) (l + 7)(2-7) 



(3.28) 



which then acquires the w-dependent shift, as explained previously. Note that Xi is the same 
as that in [11] in the n/ = limit, and differs from it by the n/-dependent terms in Agg, ^qgfl 
and Aqg. Note also that cubic and quadratic poles cancel out in xi, because of the tu-shift of 
the collinear poles and of their factorization, which are embodied in our formalism (cf. sec l4.3p . 

The remaining single poles of K^^^i are provided by which is obtained as follows. Let 
us introduce the constant (in 7) coefficients of the characteristic functions: 



Xo(7, ^) = - + Co(a;) + (7) , Xc(7, u^) = - + Uuj) + (7) . 
7 7 



(3.29) 



Then, according to eqs. fl2.22|2.23p . we have 



cAuj] 



1 



, c,(0) = 



Co(0)=0, c[,(0) 



TT 



y(i) = 



(3.30a) 
(3.30b) 



Note that Co(0) vanishes by virtue of the LLx expansion Xo{l) = 1/7 + (7^), in other words 
it is a scheme-independent coefficient. On the other hand, all other quantities in eq. fl3.30p are 
scheme-dependent, i.e., they depend on the particular choice we adopted for shifting the poles 
of xo and on the definition of Xc as in eqs. fl2.22[ 12.220 . 

By using the above notation we derive the expansions of the BFKL eigenvalue function 



BFKL 



(7) 



A 



gg 



+ 



qgfi 



72 \2Ca Ca2C, 



1 (46C^ - 52CF)Tf 



7 



9{2Ca] 



+ (7°) 



(3.31) 



the expansion of the subtraction term in eq. fl3.14p 



1 fA 



F J- gg,0 



+ 



Y \^Ca Ca 2CaJ 7 
and the expression of the pole term in Ki 

{AQCa - 52CF)Tf 



Cf) ~l" Cc 



A 



gg _^ Cp '^^qgfi 



qg 



~(0) 

Xi 



9(2Ca)^ 



+ Cc 



2C A Ca 2C a 



11 ^ {AC A - l6CF)Tf 



Ca '^Ca 



+ 0(7°) (3.32) 
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3(2Cy 



Cp S, 



qg 



Ca '^Ca 



(3.33) 



(0) 



where all quantities in the three formulas above are evaluated at = 0. Note again that xi 
satisfies the consistency relation in eq. (13.161) . proved in sec. 14.31 

The remaining expressions used in the previous subsections are easily obtained by the fol- 
lowing detailed formulas 



i^qg^gq + Tgg) + AggTggXht 



Cc^qg (Xqq + ^ gg) + ^qg^ggXht 



gq 



[K 



(i)k(o) 



igg 



CqT 



gg 



qq 

■ 2Ca 
Co 

UJ 

2Ca 



qqA 



qgA 



Cf 



+ Ce (Fgg + A 



gg) 



+ Cc \Ag(^ gg + V g^ 



gg^ 



~{0) 
^2 Al 

ijj 



(3.34a) 
(3.34b) 
(3.34c) 

(3.34d) 
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where Xht stands for Xhtil ~ 0) ~ means asymptotic in the uj ^ hmit. From the above 
formulas we can compute the high-energy hmit of Ki 



The introduction of running couphng may change our expectations on how momentum 
conservation is satisfied by our kernel. Since we incorporate NLO anomalous dimensions, as 
given by eqs. (13.221) . we expect violations at next order, that is at relative order a^. This 
is also the order at which running coupling effects start to matter in the actual derivation 
of anomalous dimensions, being related to a commutator of two values of as, evaluated at 
different scales (cf. app. Therefore, no problems arise at NLO. Since we do not explicitly 
consider incorporating NNLO results in this paper, we do not make an effort to improve energy- 
momentum conservation at frozen Og. We should add that the consistency relations on NLx 
terms which restrict our scheme start imposing a coUinear scheme restriction at order cxl/ui'^ 
which is derived in sec. 14.21 and is violated, even if marginally, in the MS-scheme. 

4 Frozen coupling anomalous dimensions 

In this section we want to discuss a number of issues concerning the anomalous dimension 
matrix in the case of frozen coupling, in which the whole matrix can be analytically calculated 
in terms of the kernel matrix elements in (7, C(j)-space. This allows us to compute the two 
eigenvalues 7 = 7±(as,ti;) and their inverses, the effective eigenvalue functions uj = x±{(^s,l), 
as well as their eigenvectors. We obtain in this way the hard Pomeron exponent co's(c^s) and the 
resummation formulae for the matrix elements of the anomalous dimension matrix. The latter, 
at a given level of the LLx hierarchy, must be consistent with the exact low order anomalous 
dimensions we have used in constructing the kernel, thus providing consistency relations for 
the collinear and fc-factorization schemes. Here we find, at NLx level, that such relations are 
identically satisfied by our construction at NLO, while they put a nontrivial constraint on the 
al/uj^ term of at NNLO. 

4.1 Anomalous dimension matrix 

If as is frozen, the matrix kernel K(A;, k'; u) is scale invarant and its resolvent admits the 
7- represent at ion 



where the ag- and a;-dependences of all the above quantities are understood. One can then 
write the spectral decomposition 





(4.1) 



where K(as,7,ct;) is the characteristic function matrix of K{k,k';uj). 

We introduce the eigenvalues ri± and eigenvectors u± of K in the usual way 



K(7)u±(7) = ri±{-f)u±{j) , 



(4.2) 



K(7) = ^+(7)n+(7)+r7_(7)n_(7), 



(4.3) 
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where Il± are the orthogonal projectors on the eigenspaces of K, and are given by 

n, = ^, n_ = ^. (4,4) 

■^±(7) being the left-eigenvectors of K satisfying = = {v-U+) when m's and v^s are 

evaluated at the same value of 7. 

The behaviour in t-space of S is determined by the 7-poles in eq. (14.11) . These poles are 
found at 7 = 7±(as,co') such that 

r7+(7+, as, uj) = 1 , r]^{-f_, a^, uj) = 1 (4.5) 

and are interpreted as anomalous dimension eigenvalues of S- 

By applying a driving term /q to the Green's function S, at leading-twist level — i.e., taking 
into account only the two rightmost poles in the half-plane 3fJ(7) < 1/2 obeying limQ,^^o7± ~^ 
— one obtains the vector of (integrated) quark and gluon densities 

(i')^f= E V^fSVnK70/o(70 (4.6) 

We want to show that / satisfies the DGLAP-type evolution equation 

d/ 



d log k? 



r/ , (4.7) 



in terms of a well-defined resummed anomalous dimension matrix F. In fact, by inserting the 
expression (14. 6 p into both sides of eq. (14.71) . the equality is satisfied provided 

[r - 7j]nK70/o(70 = , (^ = +, -) . (4.8) 

It might seem that F is dependent on the initial condition /q. This is not the case, because 
whatever the choice of /o, the projector Hiiji) projects /o into a vector proportional to ui{'-yi), 
and eq. (14. 8 p reduces to 

Tuii-fi) = -fiui{-fi) , (/ = +,-), (4.9) 

i.e., r is the (unique) matrix whose eigenvectors are {^+(7+), n_(7_)} relative to the eigenvalues 
{7+,7-}- 

P ^ _^^ -»+(7+) ^^+(7-) ^ M_(7_) ®t>_(7+) 
^ v+il-) ■ M+(7+) ~ ^-(7+) ■ W-(7-) 
In more detail, the eigenvalues 7±(q;s,u;) are provided by 

det[l-K(7±,a„a;)] =0 (4.11) 

and the eigenvectors are 

^-(7+) = (1 -P) , (4.14) 
t;+(7-) = (r 1) . (4.15) 
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Therefore, the full expression of the anomalous dimension matrix is 



1 f rp-i+ + 7- (7+ - 7-)p 



1 + ''^P V (7+ - 7-)r 7+ + rp7_ 
We obtain the relationships 

^+ = Vgg+rV,g, (4.17) 

7- = r,,-rr,,, (4.18) 

r99 = Kr,,-7-)- (4.19) 

4.2 Resummation formulae 

All the above formulas are exact in the frozen coupling case, and do not depend on our particular 
assumptions on K. Now, by taking into account the structure of K described in the previous 
sections, we compute the anomalous dimension matrix elements at NLx level. 
To this purpose, we note that the eigenvalues of K are defined by 

r72_r/trK + detK = (4.20) 

and that both tr K and det K are of order 1/a;, so that up to NLx level we have 

^+ - K - — ~ - (oKgg,o + aoKg9,i) + a iKgg^ + ^ , (4.21) 

,_:.jg^:.aY.K,.,„- °'^-";''''"° ) ■ (4.22) 



We note that the equation 77+ = 1 reduces to the usual BFKL determination of 7+ because 
o^gg,o can be replaced by a;/« in the NLx term. Furthermore, the equation r^. = 1 is dominated 
by its 7-pole part, yielding 



a 

?7_ ~ — 

7 



[iTgg,o - ^ ir,,,oJ - 1 . (4.23) 



This provides the lowest order determination of 

7_ = ^7-,o + NNLx = --^r,,,o(0) + NNLx , (4.24) 
Zn Ga 

so that, up to NLx level only the one-loop term of Tqg contributes and no small-x enhancements 
are present. 

The coefficient r is now calculable from eq. fl4.13p and, up to NLx level, we obtain 



Ca 



1 + — ri + uri + 
Ztt 



1^39,1 ~ _ Agg figg^^^Lggfi 



2Cp A 2Cp 2Ca 



(4.25) 

Note that, since eq.f l4.13p is evaluated at 7_, r does not contain l/u enhancements and is 
generally calculable from fixed order perturbation theory. The coefficient p, on the other hand, 
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is calculated at 7+ (cf. eq. f l4.12p ). and so contains resummation of NLx terms to all orders in 

Given that 7_ and p are NLx quantities, it follows that also F^g and Tgg are NLx. From the 
previous equations we obtain the resummation formulae 



F 



LLx 
91 



qg 



F 



<?<? 



^pLLx 
Ca 



(4.26) 
(4.27) 
(4.28) 



which are well known [8]. In addition, the matrix kernel predicts 



91 



Cf 

Ca l 
Cf 

Ca 



99 



11 

Cp^ 
_i 1 f_p 

ZTT Ua 



(0) + ^—n + a^n 



+ • 



(4.29) 



Note that NLx running coupling contributions are shown in Appendix A to start at order a\. 
The above resummation formula for F^^ is easily checked to be identically valid at (ctg) and 
(Og). At (ftg) it yields the relation 



pNLa; 



pNLx 



(4.30) 



which characterises the class of schemes described by our matrix formulation, and appears to 
be not satisfied in the MS scheme [26] 0, even though the violation, of relative order rif/N'^^ 
is numerically less than 0.5% for ri/ < 6. Strictly speaking, this implies that the MS scheme 
at NNLO cannot be incorporated in the present matrix approach. However, one could think of 
adding the small violation just mentioned by a matching procedure. 



4.3 Consistency relations 

They arise in general because of the joined requirements of simple cu-poles and 7-poles inposed 
on our kernel. For instance, by the 7-pole hypothesis we determine the 7-pole parts of Ki and 
K2 by the equations 

Kf) = F, - K«Fo , Kf = F, - K«Fo - K«Fi - K^^^F^ . (4.31) 

These expressions should be consistent with the w-pole hypothesis so that higher order poles 
in UJ, possibly occurring in the r.h.s. of eq. fl4.3ip . should cancel out. 

Furthermore, by the w-pole hypotesis, we determine the cu-pole part of K^^j by a subtraction 
of the BFKL kernel, as follows: 

(oKi)^^ = (2C^)2irr^L - (oKo iKo)^^ . (4.32) 

^^From eq. (4.29) of ref. [26], by taking into account the difference between our and tlieir normalization 
r^L| = -Ef/Suj^, it turns out that, in the MS-scheme, T^l^ = [r^^^ - ^] . 
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Once again, this should be consistent with the 7-pole hypothesis, so that quadratic (and possibly 
cubic) 7-poles in Kf^^^ should cancel out on the r.h.s., and furthermore the simple pole should 
be consistent with eq. fl4.3ip . that is 

(oKf))^^ = (2C^)^irr^(°) - (oKoiKo)^°; = (or,),, - oK'r,)^^ . (4.33) 

Let us start proving the consistency relation for eq. ( 14.3ip . Generally speaking, they are 
equivalent to recursive relations on the tu-singuarities of r„ or, in other words, to the resum- 



mation formulas proved in sec. 14.21 For instance, the assumed absence of a;-poles in (Kn ) 
implies the NLx resummation formulas for (r„) ^: 



qa 



(ri),.-(Ki^^),,o(ro)^,, (at level) (4.34) 



(r2),.-(Ki^\,o(ra^^, (at ^ level) (4.35) 



as predicted by eqs. KTii 14281) . 



The gq, gg entries are slightly more complicated. At order a'^/uj, oKg°\ is determined by 



eq. fl4.3ip . so that no consistency condition arises in the gq entry. However, o^^qni is already 



99, i 



determined by eq. fl4.32p . so that the consistency condition fl4.33p arises. The latter is verified 
because the simple-pole part of Kf^^^ is simply 7+1 [7], so that eq. (14.330 reduces to the 
identity (13.161) . which implies, at order a'^/u, 

7+,l - J- 33,1 - (of^O 1*^0 +0l^0 if^O )gg~\0'^0 1 '^0 +1*^0 QI^O ) gg 



oCC] =^(,K<")„, (4.36) 

Jgg U w 



where the r.h.s. reduces, by eq. (14.340 . to ^r^^ i, as it should. 

Furthermore, at order ag/a;^ we have consistency conditions for Tgg^2 and Tgq^2- The former 
is identically satisfied, by some algebra similar to eq. (14.361) . because 7+1^"^ = ^^^,2 + ^^qg',2 
given in eq. (I4.17p . The latter is instead non trivial and, after a similar algebra, reduces to 



J- 99,2 - ^ J- 33,2 ) l'i-<3'J 



as already proved in eq. (I4.30p . with the same consequences. 

We finally note that cubic and quadratic 7-poles are absent in (I4.32p because of the identity, 
valid up to order 1/7^, 

-(«Kf)„(,Kf)„+(„C)^,(,Kr),,. (4.38) 

Here the cubic poles at 7 = already cancel out in the l.h.s., because of the |xoXo subtraction 
needed to switch energy-scale kko —>■ fc^, due to the a;-shift (12. 8p . The remaining quadratic 
poles are given by the r.h.s., because of normal coUinear factorization, and of absence of I/7 
singularities in xq at energy-scale k"^ (eq. (I3.15P ). 
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5 Characteristic features of the resummed Green's func- 
tion 



In this section we present numerical results of some plienomenologically relevant quantities 
which can be obtained by using the matrix kernel K = dKo + a^Ki developed in the previous 
sections. We recall that the final form of Kg and Ki can be found in eqs. (I3.2p and ( 13.220 
respectively, and that the detailed implementation of the running coupling is found in eq. fl3.27p . 

We state once more that our matrix kernel incorporates exactly the DGLAP and BFKL 
properties at NLO and NLx accuracy. However, in order to see the impact of the NLO contri- 
butions and to compare with previous resummation approaches, we will consider also results 
obtained from the kernel with only LO anomalous dimensions (but still in NLx approximation). 
The corresponding kernel — which we refer to as NLx-LO model — is built with the same Kq 
given in eq. (13. 2p but with Ki including only the Ki term in the gg entry, as can be read from 
eq. fl318D by setting Ti = 0. 



5.1 Hard Pomeron exponent 

We shall first investigate the high-energy s — +oo behaviour of the A + B ^ X differential 
cross section given in eq. (12.61) at fixed and equal value of the two hard scales ~ k"^, by 
determining the growth exponent (hard Pomeron) a;s in the limit of frozen coupling. In this 
limit, we can use the representation (14. ip for the Green's function S{k, k'; u) and, by using the 
spectral decomposition introduced in sec. 14. we obtain (the ag-dependence is understood) 

daik^ko) ^ /• d7 f duj f s y h^{k,uj)Ui{j,u)h''{ko,u) 

- . (5.1) 



The a;-integral gets contributions from the singularities (labelled by the index m) of the inte- 
grand at = uji^mil) due to the vanishing of the denominator 



r]i{^,ui,U-f))=0, (5.2) 



thus providing 



da(A;~A;o) /" ^7 / s V'-^^^ /i^(fc, cD) (7, c^) /i^(fco. 



d^kd^ko ,^f^^_^^J 27ri \kkj -dM^,u) ■ ^ ■ ' 

In the limit s ^ kko the 7-integral is dominated by the saddle-point 7 = 7s such that 

-^(^i,mhs) = (5.4) 
d7 

for the particular values of / and m such that uii^mi'^s) is maximum. It turns out that those 
values correspond to the leading-twist component of the (/ = +)-branch of the eigenvalue 
function rj^, namely the solutions of ujs) = 1 with 7s — * for — ^> 0. As a result, in the 

high-energy limit the cross section has the power-like behaviour 

dajkc^ko) _ ( s Y 
d^kd^ko ~ \kkj ^^-^^ 
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where the process dependent coefficient Cg is constant or at most logarithmic in s, and the 
growth exponent a;s is determined by the conditions fl5.2ll5.4l) . 

We can recast eq. fl5.4p into an equivalent relation for the function d^rji. In fact, by taking 
the total 7-derivative of eq. fl5.2p we can express 

thus obtaining the following conditions for the hard Pomeron exponent: 

r]+{%,uJs) = l (5.7a) 
djr]+{%,ujs) = . (5.7b) 

The above conditions in turn can be translated into analogous conditions for the determinant 
of the operator 1 — K. In fact, from the relations 

det(l - K) = (1 -r/+)(l -r7_) (5.8a) 
d^det{l - K) = -[{1 - r]+)d^r]_ + (1 - r/_)9^r/+] (5.8b) 



eqs. (15.71) are equivalent to 



det[l - K(7s,cui)] = (5.9a) 
a^det[l - K(7s,cui)] = . (5.9b) 



0.5 




Figure 2: Hard Pomeron exponent obtained in the NLx-NLO matrix formulation with Uf = A 
(solid blue) and Uf = (dotted blue). The one-channel results [11] are also shown (red squares) 
and compared to those of the matrix model in NLx-LO approximation with nj = (solid green). 
The calculation is done in the fixed coupling case. 

We have numerically solved the implicit equations (15.91) in our matrix formulation in a range 
of as up to 0.4, both at NLx-NLO and NLx-LO accuracy, for two values of n/ = and 4. The 
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results for ujg versus as are shown in fig. [2|, where we compare with results obtained from our 
previous one-channel approach. 

The NLx-LO curve at rij = almost overlaps to the old one-channel result, thus showing the 
stability of the matrix formulation and the continuity with the one-channel approach. In fact, 
since at rij = the kernel is diagonal, only the gg entry determines tUg. The small discrepancy 
is due to: (i) the momentum-conserving factor in front of Xi in sq. (13.181) (it was a plain l/u 
in the one-channel case); (ii) a non vanishing two-loop anomalous dimension (only for the low- 
energy part, actually) (we enforced vanishing anomalous dimension in the one-channel case). 
By including the NLO contributions we obtain a moderate increase of the Pomeron intercept, 
which slightly diminishes when quarks are also taken into account. 

5.2 Effective characteristic function(s) 

As already noted in the previous section, the contributions to the integral representation of the 
cross section stem from those values of u and 7 such that 



which provides a relation between the moment index 00 and the anomalous dimension variable 
7. Solving eq. (15.101) for either or 7 defines the effective characteristic function and its dual 
effective anomalous dimension 



While in the one-channel case we have only one perturbative branch of those functions, cor- 
responding to the BFKL eigenvalue function x+{l) ^ind to the larger eigenvalue 7+(a;) of the 
anomalous dimension matrix, in the matrix formulation we expect two branches. The sec- 
ond branch corresponds to the smaller eigenvalue 7-(c<j) which is dual to a second effective 
characteristic function X-{l)- 

In fig. [3] we show the two branches of the effective characteristic functions obtained in the 
NLx-NLO and NLx-LO cases. We have considered here the asymmetric cj-shift corresponding 
to the energy-scale fc^, with as = 0.2. The x+'s are characterised by the typical minimum 
around 7 ^ 0.5 whose value is nothing but ti;s(as). On the other hand, the X-'s appear as 
steeply decreasing functions located around 7 < in the region shown in our plots. 

The continuity of the resummation procedure when going from the one-channel to the two- 
channel formulation at = is illustrated in fig. [3^ by the overlapping of the (-l-)-branch of 
the NLx-LO curve to the circles corresponding to the one-channel scheme-B effective eigenvalue 
function. The NLO terms provide a slight increase of the (+)-branch in the region < 7 < 2, 
and a small decrease of the (— )-branch at 7 < 0. At = there is a crossing point of the two 
branches at negative 7 in either approximations. 

The inclusion of quarks removes the crossing (with a mechanism similar to the degenerate 
level splitting in quantum mechanics) causing X- to be always on the left of x+5 as can be 
seen in fig. [3)3. Quantitatively, the quark contribution lowers both x+ in the region around the 
minimum (compare the two inserts in fig. [3]) and X-- 

Note the two fixed points at (7,0; = Xeff) = (0, 1) and (2, 1) of the (+)-branches. In the 
one-channel case these fixed points corresponds to momentum conservation in the coUinear and 
anti-coUinear limits respectively. In the two-channel formulation they imply that the anomalous 
dimension eigenvalue 7+(co' = 1) = 0; however, momentum sum rule is satisfied provided the 




(5.10) 





(5.11) 
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Figure 3: Effective eigenvalue function obtained in the matrix formulation with full NLx-NLO 
accuracy (blue) and with NLx-LO accuracy (green) fornj = (a) and n/ = 4 (b). The steeply 
decreasing curves on the left side of each plot represent the minus-branches X-, while the curves 
with a minimum around 7 ~ 0.5 represent the plus-branches x+- The red circles reproduce the 
one-channel result in scheme B of ref. [11]. The calculation is done at fixed coupling = 0.2. 



as = 0.2 nf = 



NLx-LO 
NLx-NLO 
1 -channel B 
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corresponding left-eigenvector v+ = (^r l) of the anomalous dimension matrix (cf. eq. (14.141) ) 
be (l l) at cij = 1, i.e., provided r{uj = 1) = 1. 

Actually, our matrix model presents a small violation of the momentum sum rule. In fact, 
by exploiting the fact that 7_|_(u; = 1) =0 and by using eqs. (I2.19p and (14.161) we have (at 

UJ = 1) 

q + g=il l)rf^^ =7-^^(g-p^7). (5.12) 



1 + rp 

The computation of the prefactor 7_ (1 — r)p/{l + rp) versus ctg shown in tab. [1] gives us an 
estimate of the relative amount of momentum non-conservation; the violation is of order for 
the NLx-LO scheme, and of order for the scheme with NLO terms included. 





NLx-LO 


NLx-NLO 


NLx-LO/a^ 


NLx-NLO/a^ 


0.025 


0.00019 


0.0000031 


0.302 


0.199 


0.050 


0.00072 


0.0000208 


0.287 


0.167 


0.100 


0.00260 


0.0001303 


0.260 


0.130 


0.150 


0.00534 


0.0003437 


0.237 


0.102 


0.200 


0.00872 


0.0006107 


0.218 


0.076 



Table 1: Estimate of momentum sum rule violation. The quantity '-f^ (1 — r)p/(l + rp) at uj = 1 
has been computed for various values of as (column 1) in the NLx-LO (column 2) and NLx- 
NLO (column 3) schemes. Column 4 (resp. 5) shows that the NLx-LO (NLx-NLO) violation 
is of order a"^ (c^D- 



6 Numerical results with running coupling 

In this section we shall present results obtained by solving eq. (12. ip in {x, fc)-space, including a 
running coupling. The basic structure of the ensuing integral equation follows from eq. (13.270 
and reads (Y = logl/x) 



6{e - kl) 



7r 



+ E / — / a{q^)6ag6gcK,{z-k,k') 



+ «(fc2 )Kcoii,ac(^; k, k') + a\kl)K^,ac{z- k, k') gefe(log -; k\ ko) (6.1) 

(the explicit expressions of the kernels Kq, Kcou and Ki in the equation above are given in 
app. E]). We shall extract Green functions and sphtting functions, using the methods described 
in [11, 12,23]. In both the coupling and the kernels we use a fixed number of flavours, Uf = 4. 
The coupling runs with a 2-loop f3 function, and is normalised such that as{3 GeV) = 0.256. The 
infrared region of the coupling is regularized by setting it to zero for scales p < po = 0.75 GeV. 

The results that we shall show are those of the model described above (NLx-NLO), and 
also those for a model in which the higher twist part of K^^gg has been supplemented with 
(sjTiimetric) 1/(1 + 7)^ and 1/(2 + u; — 7)^ terms [£| so that not only the al/u but also a^/u'^ 
terms of the qq, qg and gg splitting functions are in the MS scheme!^ We shall denote this 



i^More precisely, the higher-twist kernel reads XhtM = ^(1 + 7)"^ - if (1 + 7)"^ + [(1 + 7) ^ (2 + t^-7)]- 
^^The gq term is almost in the MS scheme, the only difference being a small A^c-suppresscd contribution of 
relative order Uf/N^, corresponding to the violation of eq. (|4.37p in the NNLO MS splitting functions. 
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Figure 4: Gluon-induced part of the Green function for the NLx-NLO and NLx-NLO'^ models, 
compared to the results of [12] (scheme B). For the models of this paper both Sgg and Sqg are 
shown. The value chosen for the coupling, as = 0.15, corresponds to ko ~ 20 GeV. The band 
indicates the spread in the result for the NLx-NLO model when varying the renormalisation 
scale in the range 0.5 < x^ < 2. 



second model NLx-NLO"*". We shall also compare to results obtained in our earlier single- 
channel work [12] (scheme B), where we used a 1-loop coupling, nj = in the kernel (but 
n/ = 4 in the coupling) and for which the NLO piece of the effective Pgg splitting function was 
identically zero. 



6.1 Green functions 

The Green function for the matrix evolution is itself a matrix in flavour space. Physically the 
most interesting part is that involving gluonic sources, and this is shown in fig. |H 

The old and new resummations give nearly identical results for the Sgg part of the result, 
indicative of the stability of the resummation procedure. Furthermore, the differences between 
them are much smaller than renormalisation scale uncertainty, which grows with Y. The 
growth with Y of the scale uncertainty can be understood as an indication of underlying scale 
dependence of the effective BFKL exponent. 

The Sqg channel can be given only within the new resummations. As would be expected, 
the quark component is suppressed by a factor ~ as compared to the gluon component. A 
consequence of the fact that the quarks are only generated radiatively is the scale dependence 
in their normalisation as well as in their growth with Y. We note that the change induced 
by the NNLO scheme-dependent higher-twist part of Ko^qg (NLx-NLO"*" versus NLx-NLO) is 
small, despite the fact that the region A; ~ /cq that we study is that most likely to be sensitive 
to this higher-twist contribution. 
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Figure 5: The matrix of NLx-NLO (and NLx-NLO^ ) splitting functions together with their 
scale uncertainty and the NLO splitting functions for comparison. In the gg channel, we also 
show the old scheme B result (rif = 0, no NLO contributions, 1-loop coupling) of [11]. The hand 
corresponds to the span of results (NLx-NLO) obtained if one chooses = 0.5 and = 2.0. 



6.2 Splitting functions 

The extraction of splitting functions is carried much in the same way as in the one-channel 
case described in [11,23]. There a special (infrared) inhomogeneous term was included in 
the equation for the Green function such as to ensure that the resulting integrated gluon 
distribution satisfies ^(^(x, /i^) = 1, independently of x, for /i^ set equal to some given Q^. 
With that inhomogeneous term fixed, the xPgg{x,Q'^) splitting function was then obtained as 
^^■g^-^xg{x, ^'^)\^2^Q2. In the matrix case, we have a 2-component vector of inhomogeneous 
terms: we can choose it such that = 0, xg{x,fi'^) = 1 for /x^ = Q^, in which case we 

obtain 



xPgg \ d d 



xq{x, jj?) 

2^ 



xPgg J 9 In a; (9 In /i2 \xg{x,^i^) J ^^=^2 ' ^^'^^ 

Alternatively we can set the inhomogeneous terms so as to ensure that q{x, yU^) = 1, g{x, /i^) = 
for /i^ = and we then extract Pqg and Pgg. 

The matrix of effective splitting functions as determined with this method is shown in fig. [5l 
for both our kernels and with a scale Q ~ 6 GeV, giving as(Q^) = 0.2. For reference we plot 
also the exact NLO splitting functions and our previous results for the single-channel evolution. 
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Considering first the gg channel, the results are rather similar to the old ones, and in particular 
maintain the characteristic dip [13] around x = 10~^ that has been seen also by the authors 
of [16]. This dip is present also in the gq channel, and indeed the gq channel is rather similar 
in a range of features to the gg channel, which is natural since it is largely driven by the 
summation of the g g branching. Among the common features is the slight but noticeable 
difference compared to the NLO splitting functions at moderate x (x > 0.1). The detailed 
origin of this characteristic is not really understood, but may well be connected with the fact 
that the various pieces of the NLO gg splitting function are effectively placed in different parts 
of our evolution kernel, and then subjected to non-trivial (higher order) non-linear effects in 
their recombination into the final effective splitting functiono 

Concerning the scale dependence of the Pgj channels, as for the Green function, it grows 
significantly towards small x, and again this is a sign of scale dependence of the small-a; in- 
tercept. One notes non- negligible scale dependence also at moderate x. Generally, down to 
moderately small x, the scale dependence of the NLx-NLO splitting functions (all channels) is 
rather similar to that (not shown) of the plain NLO splitting functions. 

The two Pqj channels differ fundamentally from the Pgj channels in that they are non-zero 
at small x starting only at NLx and at NLO. Thus there is a sense in which our NLx-NLO 
treatment is effectively a leading order treatment for these channels, at least as concerns their 
normalisation (the small-x growth is driven by iterations in the gluon channel, so one expects 
this to be under better control). This is visible in the much larger scale dependence for these 
channels. They also have some (modest) sensitivity to the difference between the NLx-NLO 
and NLx-NLO^ kernels, whereas in the Pgj channels there was almost no sensitivity to this 
difference (even though the difference is NLx in all channels). A general feature of the Pgj 
splitting functions is that they are rather similar to the NLO splitting functions (more so than 
in the gluon channel). In particular, though like the Pgj splitting functions they have a dip 
around x = 10~^, this dip is considerably shallower. The conclusion here is that the NLO 
Pqj splitting functions can probably be considered a good approximation to the full splitting 
functions for x as low as 10~^. 

An important cross-check of the methods used to extract the splitting functions is that the 
results should be independent (modulo higher-twist contributions) of the infrared regularisation 
of the coupling, i.e. independent of the scale fio below which the coupling is set to zero. To 
this end we have extracted the splitting functions with /xq increased from 0.75 GeV to 1 GeV 
(corresponding to reducing as(/io) from 0.58 to 0.45) and find that the results change only by 
a few percento As in previous work [11,23] we find that these factorization violations scale 
roughly as 1/Q rather than as l/Q^, a characteristic perhaps attributable to resummation 
effects, which could quite conceivably modify typical collinear l/Q"^ power-suppressed effects 
such that they become 1/Q'^~'^'^ with an effective 7 ~ 1/2. 

We close this section by showing in table [2] the degree of momentum sum-rule (MSR) 
violation in the splitting functions for three values of as- From just a small number of values 
it is difficult to establish the exact scaling law0 and in particular it is difficult to determine 
the relative admixture of higher-twist and perturbative components in the MSR violations. 

^^The greater similarity between the large-x scheme B kernel and the NLO results is an artefact related to 
the different Uf values used in the old scheme B results and the new matrix evolution. 

"'^^One may also reduce //q, however for /xq ^ 0.6 GeV, as(/xo) then becomes so large that numerical instabilities 
develop, and it becomes impossible to extract meaningful results. 

^^Limits on the available numerical accuracy make it difficult to obtain reliable estimates of the MSR violations 
for smaller values of as, because as as decreases one needs ever higher relative accuracy to accurately determine 
the rapidly vanishing MSR-violating component. 
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NLx-NLO 


NLX-NLO+ 




Q [GeV] 


E,r,,(i) 


E,r,,(i) 


E,r,,(i) 


E,r,,(i) 


0.20 


6 


0.0079 


-0.0059 


0.0074 


-0.0055 


0.15 


20 


0.0021 


-0.0015 


0.0018 


-0.0012 


0.10 


220 


0.00012 


-0.00003 


0.00006 


0.00002 



Table 2: Momentum sum-rule violation in the NLx-NLO and NLx-NLO^ models for three 
values of as- The numerical uncertainty is roughly ±1 on the last digit of each result. 

Nevertheless, one sees that the MSR violation vanishes very rapidly as ag decrease, suggesting 
that a significant component of it is non-perturbative in origin. This conclusion is borne out 
by studies which show that the amount of MSR violation depends somewhat also on /ig, the 
infrared cutoff scale for the coupling. 

7 Discussion 

We have proposed here a matrix evolution equation for the flavour singlet, unintegrated quark 
and gluon densities, which generalizes the DGLAP and BFKL equations in the relevant limits. 

The matrix approach (sees. [2] and [3]) is supposed to unify collinear and high-energy factor- 
izations in both partonic channels, and is not necessarily guaranteed to actually work, because 
of the various crossed consequences that the above factorizations have: consider, for instance, 
the anomalous dimension resummation formulae arising from fc-factorization [5, 8] and the 
7 ^ 1 + a; — 7 symmetry of the BFKL kernel [9, 10] arising from collinear factorization. It 
is therefore a nontrivial result of this paper that our resummed splitting functions do satisfy 
collinear factorization in matrix form, as shown in sees. |l]and[6l In this respect, our approach 
defines, by the matrix evolution, some unintegrated densities that are appropriate both in the 
collinear and in the small-x limits. It would be interesting to explore the relationship of such 
explicit construction with alternative studies [22,24]. 

Furthermore, we want to incorporate exact low-order anomalous dimensions in our matrix 
kernel, say in the MS scheme. We find, in this context, a new kind of consistency relations on 
the kernels, due to a possible clash of exact low-order expressions with a novel NLx resummation 
formula for F^^, arising in the matrix evolution (sec. Hj). We prove such relations to be satisfied 
by our construction in the MS scheme at NLO, but marginally violated by ny/A^'^-suppressed 
terms at NNLO. We are thus able to complete our construction with exact NLO anomalous 
dimensions and NLx kernel, and we postpone the analysis of the NNLO accuracy, which is 
however nearly incorporated (in the NLx approximation) in our NLO"*" version. 

The frozen-Os features of our matrix model are characterized by the previously mentioned 
resummation formulae of sec. HI and by the hard Pomeron exponent and effective eigenvalue 
functions of sec. [5l One should notice the basic continuity of our matrix approach with the 
single-channel case in the = limit, and the corresponding agreement of the leading effective 
eigenvalue function with the ABF approach. Additionally, we provide here the subleading 
effective eigenvalue at = 4, corresponding to the 7_ eigenvalue of the anomalous dimension. 

We are finally able to provide the whole matrix of resummed splitting functions in sec. El 
Roughly speaking, the outcome shows that resummation effects are small in the Pqa entries up 
to a;-values as small as 10~^, while the shallow dip is the main qualitative feature of both Pga 
entries, with resummation effects starting below x ~ 10~^. 
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We notice that the above results are in the MS scheme up to NLO (including approximately 
NNLO in their NLO"*" version), but generally differ from it at higher orders. However, resum- 
mation effects in the scheme change have been studied in [14] and turn out to be of the order 
of the scale uncertainty. For this reason we believe that the above results can be safely used 
in the study of structure functions and other cross-sections, by supplementing them with the 
corresponding coefficient functions or impact factors. 
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A Recursive expressions for the anomalous dimensions 

In this appendix we show how eqs. fl2.20p are obtained. Let us first rewrite eqs. (12.1 8112. 19p in 
7-space, by noting that 7 = 5iogQ2: 

7 Ml, ^) = ^iil, ^) = ^ij{^)fj{l, ^) + h.t. (A.l) 

where "h.t." stands for higher-twist contributions characterised by being regular at 7 = 0. It 
follows that, in matrix notation, 

7'/ = 7(r/ + h.t.) = r(7/) + h.t. = + h-t. (a.2) 

and, by induction, 

7"/ = r"/ + h.t. , (A.3) 

Secondly, we consider eq. (I2.17P and expand the matrix kernel K in powers of gamma (according 
to the notations following eq. (12. lip ), obtaining 

^= [Yl K('")7'"-i j + J^"""""" = k('")7'"/ + jrsource ^ ^ K^^^^P'"/ + h.t. , (A.4) 

\m,=0 / m=0 m=0 



By comparing eqs. flA.ip and flA.4p we derive the implicit equation 

r = Y K('")P'" , (A.5) 

m=0 

which allows us to determine the effective anomalous dimension matrix P in terms of the matrix 
kernel K. 

It is now straightforward to compute the perturbative coefficients P„ defined in eq. (I2.12p . 
By expanding eq. flA.SP to first order in a yields (m = 0) 

'(0) 



dPo = aK^o' (A.6) 
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from which eq. fl2.20al) follows. By expanding eq. flA.51) to second order in a yields (m < 1) 

dFo + d^Fi = oKf + 6?Y.f + dK^^^dro . (A.7) 

At frozen coupling, the operators Kq^^ and a commute. By then collecting the (d^) terms 
and remembering that Fq = Kq'^^ we get 

c?V^ = a\Y.f ^Y.^^'^Y.f) (A.8) 

from which eq. (]2.20bp follows. A similar iteration procedure produces eq. fl2.20cp and higher 
orders. 

In the running coupling case, we have an additional commutator term starting at second 
order, namely 

d[K[,'\d]K[,°^ (A.9) 

By using the expansion 

a = d^- Po&l log g + (dy , (A.IO) 

the commutator reads 

[K(^\d] = -/3od2[K«,logg] + 0(d3) , (A.ll) 

thus producing, by eq. flA.Qp . a contribution of order d^ to the anomalous dimension matrix. 
Note however that the NLx term of order ag/cu^ vanishes, because the gg entry of Kq^'* has no 
leading l/u term. 

Extending the above procedure to higher orders, we see that at each order (say N"LO) 
the anomalous dimension gets a new term Kn^ but also a series of other terms which are 
combinations of the anomalous dimensions at the lower orders < n. 



B Splitting functions and anomalous dimensions 

Singlet anomalous dimensions and splitting functions at lowest order appear in many of our 
formulas. The former are given by {Tf = TrUf) 



Pqqfl{z) — Cp 



2 I /I 



l-z + -5{l-z) 



P,g,,{z)=2Tf[z' + {l-z 

Pgqfl{z) = Cp — 

P,,,oiz) = 2Ca 



z 

1 1 



z — z 



6 



6{l-z) 



(B.l) 
(B.2) 
(B.3) 
(B.4) 



The anomalous dimensions, i.e., the Mellin transforms of the splitting functions, are given by 

1 1 3 

7TT " uj + 2 ^ 2 

1 



^qq,o{^^) — Cf 



2iIj{1) - 2iIj{uj + 1) 
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UJ U + 1 07 + 2^ 

^1 . X . X 2 1 
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oj+l ' UJ + 2 ~ cu + 3 ^ 12 



2Tf 



(B.5) 
(B.6) 
(B.7) 
(B.8) 
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The gluon anomalous dimensions Tgg^ and r^^^o are singular at u; = 0. The regular parts are 
defined by subtraction of the cj-pole, namely 



Their values at = are 



2Cf 



UJ 



^93 (^) = rg3,o(^) - 



2C, 



%,,o(o) r,,,o(o)^ 

Ag,{0) Agg{0) 



^ 

^ 3 



3Cp 



11Ca+4Tj 
6 



(B.9) 
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C Kernels and characteristic functions 

Our method of resumming energy-scale dependent terms relies on the introduction of improved 
kernels whose characteristic function^ Xcoil) are cj-dependent. In general such characteristic 
functions are symmetric in the 7 — > 1 +u; — 7 transformation and have the following structure: 

-xM = M(^) [xM + + ^ - 7)] , (c.i) 

UJ 

where the left-projection yt, contains coUinear (and possibly higher-twist) singularities only in 
the half-plane 3^(7) < 00 

The cu-dependence in the argument of the second xl term imposes kinematical constraints 
on the longitudinal momentum fraction variable z conjugated to u. In fact, by denoting by 
zM{z) and Kl the inverse Mellin transforms of M and xl respectively, we have 



"^'■'"'"'^ (V^y = .'M(.')A-.(»:>.*<) . 



(C.2) 



where fc< = mm{k, k'), ky = max(A;, A;') and z' = z ■ max(l, A;'^//c^). Since < z' < 1, the 
kinematical constraint k'"^ < k'^/z follows. 

The lowest-order matrix kernel Kq in {z, fc)-space can be derived from eq. (13. 2p and reads 



[q;Ko](2;; k, k') = a{q^) 
+ a{kl) 



^0 

,0 2CAKo{z;k,k') 



yPgqfi{z') Pggfi{z') - ^ 



Kcik, k') + z' 



(C.3) 

^0 ^,g{z')K^,{k,k') 



.0 







having defined 



Kc{k, k) — — , 
Ku{k,k')^-^, 

Agg{z) =5qg3{l 



(C.4) 

(C.5) 
(C.6) 



^^Apart from the (running) coupling factors, we always deal with scale-invariant kernels. 

^^Note that in this article we adopt the asymmetric — upper in the notations of [11] — energy scale sq — k^, 
since fc^/s is the correct Bjorken scaling variable in the coUinear limit k ^ k' we are interested to. This causes 
the ijj-shift to apply only on the (1 — 7) argument of xl and asymmetric kinematical constraints in the z variable 
as shown below. 
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The terms in square brackets in eq. flC.3l) correspond to the operator Kcou introduced in eq. (16. ip . 
The action of the first term ~ a(g^) on a test function f{x, k) is 

£ ^ 1 dk" a{q')Ko{z; k, k')f{-^, k') (C.7) 

^ t- maiq')\f{-,\k + qme-zk'')-Qik-q)f{^,k) , 

while the action of the terms ~ a(/c>), e.g. the higher-twist one, is 

f'^ I a{kl)z'^,g{z')Ku{k, k')f{-, k') (C.8) 

^1 dfc'^ &{k')z^,,{z) + j^^ dk'^ «(fc'2)^|_A^^(^|_) A:') . 

In order to obtain Ki(z; fc, k') according to eq. (13.221) . we start by computing the first term 

proportional to Xc K^.. The inverse Mellin transform of r^'^^'' is just the matrix of the 
two- loop singlet splitting functions in the MS-scheme [27]. The inverse Mellin tranform of the 
subtraction Kq^'*Kq°'* can be either computed by inverting the expressions listed in eq. (I3.34p . 
or by convolution in z-space of the corresponding factors. Here we choose the second method, 
by computing first the analytic expressions of all factors in eq. (13.340 and then the numerical 
convolution of the ensuing functions. We already obtained the Mellin transform of Ag^ in 
eq. flC.6l) : the transforms of the Tij factors are just the one- loop splitting functions reported in 
eqs. flB.5tiB.8|) : the remaining functions are listed below: 

p,z--cM = z, (C.9) 

^ Z-- ^Co(a;) = 2C^log(l - z) , (C.IO) 



J ^^^~^Xu{0,uj) = ^[6{l-z)+z'] . (C.ll) 



Finally, we need the inverse Mellin transform 



Z-- (--^]=l~2z, (C.12) 



27ri \uj 1 + UJ 
and the kernel 

K,,,, = K,-xf^K,, (C.13) 
whose a;-shifted form occurs directly in eq. (I3.27p . and has characteristic function 

XMeg = XI - xf^Xc (C.14) 
= Xi- xT - Xo ( Xo + Xc 1 - ^Xc ( Xc + Xht^ 1 - Xi Xc • 

The numerical coefficients Agg{^)^Vqg^Q{Q) can be found in eq. fIB.lOp . 5qg = Agg(O) in 
eq. (13. 8p and Xi^ in eq. (I3.33p . Furthermore, the computation of the kernel -ft"i,reg requires the 
subtraction from the NLx BFKL kernel [6, 7] of the running coupling terms and of additional 
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kernels corresponding to the characteristic functions on the r.h.s. of eq. flC.14p . They are given 
by 



Xc 



XoXo 



XoXc 



XcXht 
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logTT 



kl 
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(C.15) 
(C.16) 
(C.17) 
(C.18) 
(C.19) 



The resulting expression for -ft'i^reg is 
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(C.20) 



where (-ft'o) denotes the azimuthal average of the LLx BFKL kernel whose action on a test 
function /(fc) is given by 



&k 



12 



\k"-k^\ 



2kl 



k'^ + k^ 



fik) 



(C.21) 



Finally, we provide the eigenvalue function of -ft'i,reg with kinematical constraints, which is 
given by Xi,rcg,L(7) + Xi,reg,L(l + ^ — l), and occurs directly in eq. (]3.27p . The left projection 
Xi,reg,L of the eigenvalue function (1C.14I) can be computed starting from the expression of xil 
in eq. (A. 13) of ref. [11] and noticing that: 

(i) here we have more terms to subtract, namely those proportional to Cf/Ca and xf''; 

(ii) in ref. [11] we subtracted the n/-part of the double pole by letting Agg{0) — >• 2CaAi{0) = 
Agg{0) + {CF/CA)^ggfi{0), while here we keep only Agg{0) in front of XoXc since the n^-dependent 
double pole is subtracted by the Tggxl term. 

Therefore, in order to complete the calculation we need the left projections of the following 
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kernels: 



[xl]L{l) = [XcL{l)f + 2XcL{l) 



1 2 

— + - 

7 7 



[XcXht]L(7) = Xcl(7) - ^XhtL(7) = - - ^TT^ — r 

2 7 3(1 +7) 



(C.22) 
(C.23) 



The final result is 



Xl,reg,L(7) = [V^(l)-V^(7)] 
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(C.24) 



where 



n(7) = dtf 



.,Li2(l)-Li2(t) 
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n + 7 
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-^'Q)+^'(7)-^'(7 



(C.25) 
(C.26) 
(C.27) 



The explicit form of -f^i^reg in {k,z) space is obtained by use of eq. (1C.2I) . or by introducing 
the kinematical constraints on eq. ( ]C.20I) directly. The final form of Ki{z;k,k') follows from 
eq. dMZD- 
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